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ABSTRACT 


It is shown that the contribution of higher harmonics to 
the lunar transfer functions for the tangential components of 
the surface magnetic field is significant at frequencies greater 
than 0.01 Hz (Apollo 12-Explorer 35 magnetometer data extend to 
frequencies as high as 0.04 Hz). The inclusion of the higher 
harmonics shows that there are two distinct transfer functions 
corresponding to the components of the tangential surface 
magnetic field perpendicular and parallel to the direction of 
the wave vector of the external disturbance forcing the lunar 
induction. The dependences of these transfer functions on 
frequency and location are determined. The effects of the 
higher harmonics can 1) account for a hitherto unexplained 
feature in the Apollo 12-Explorer 35 transfer functions, namely 
the rolloff at high frequencies and 2) offer a possible explana 
tion for the frequency dependence of the difference between the 
transfer functions for the two orthogonal components of the 
surface magnetic field. The harmonic response of a simple 
current layer model of the Moon is derived and shown to provide 
a reasonable fit to the experimental data. 
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INTRODUCTION 

The Apollo 12 Lunar Surface Magnetometer (LSM) experiment 
has stimulated considerable theoretical investigation into the 
interaction between the Moon and magnetic field fluctuations in 
the solar wind. The Moon-solar wind magnetic field interaction 
can be characterized by transfer functions for the components 
of magnetic field tangent to the lunar surface, i.e. the ratios 
of the tangential magnetic field components at the lunar surface 
to the tangential components of the solar wind magnetic forcing 
field. Previous approximations to these transfer functions [Blank 
and Sill, 1969] are valid only at low frequencies and thus have 
been unable to account for the high frequency rolloff of the 
experimental lunar transfer functions (see Figur'e 1, in particu- 
lar the frequency range 0.01 - 0.04 Hz). It will be shown that 
at high frequencies the transfer functions for the two orthogonal 
tangential magnetic field components are in fact different, 
dependent on position on the surface of the Moon, and can exhibit 
rolloff at high frequency. The behavior, with frequency, of the 
difference between the transfer functions for the two orthogonal 
components of the tangential magnetic field can be understood 
by the high frequency effects to be discussed in this paper. 

In the following section we outline the theory leading to 
the definitions of the tangential magnetic field transfer functions 
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and discuss their relationship to the low frequency approxi- 
mations. Transfer functions are then presented for several 
lunar electrical conduct i v i ty model s and are qualitatively 
compared to the experimental data. 

LUNAR RESPONSE TO INDUCTION 
BY THE SOLAR WIND 

A complete solution of the electromagnetic interaction 
between the solar wind and the Moon requires that one solve a 
complex magnetohydrodynamic flow problem external to the Moon 
coupled w.ith the matching at the 1 una r' surface of the electro- 
magnetic fields in the Moon's interior. An approximate solution 

for the sunlit side of the Moon can be obtained because of the 

* 

experimental observation that the solar wind flow shows no 
marked perturbations except for the diamagnetic cavity and the 
apparently minor after body shocks at the cavity-solar wind 
boundary [Colburn et al., 1971]. To a first approximation, the 
effects of the cavity are negligible on the sunlit lunar hemi- 
sphere where, in addition, the dynamic pressure of the solar 
wind is effective in confining induced magnetic fields to the 
interior of the Moon [Sonett et al., 1971a, b] . Thus as in 
previous analyses [Blank and Sill, 1969; Schubert and Schwartz, 
1969], the electromagnetic Moon-solar wind interaction problem 
for the sunlit hemisphere will be solved by assuming that the 
Moon is immersed within a perfectly conducting space. Physically, 
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the induced dayside field will be confined by an external 
current layer above the lunar surface. By placing the 
perfectly conducting space at the lunar surface, it is assumed 
that the thickness of the confining current layer is small 
compared to the depth at which lunar induction currents flow. 

The adequacy of this approximation can be assessed by noting 
how closely the transfer function for the component of magnetic 
field normal to the lunar surface approaches un i ty- [ Sonett et al . , 
1971a]. The Apollo 15 subsatellite magnetometer, which will 
sample the magnetic field at altitudes of less than 100 km above 
the lunar surface, may be expected to clarify the nature of the 
confining current layer. 

Following the development in Sonett et al . [1971a], the 
solar wind forcing magnetic field oscillation is written in the 
form 


H = t]Hq exp 2ni(CA - ft) 


(1) 


where the Cartesian coordinate system (§,T),C) with unit vectors 

A A 

( 5 » t) » C ) 's fixed in the Moon with its origin at the Moon's 
center (Figure 2). The quantities Hg , X and f are the amplitude, 
wavelength and frequency, respectively, of the magnetic field 
oscillation propagating in the positive C - d i rect i on . 



Sonett et a] . [1971a] have shown that the lunar magnetic 

response is primarily a transverse electric one with magnetic 

field components (for the total magnetic field) which can be 

■■ 2rr i f t 

written in the form (aside from the factor e ) 
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where (r,9,cp) are spherical coordinates with the polar axis 
along £ and r=a is the lunar surface (see Figure 2), are the 
spherical Bessel functions, 
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p| (cos9) are the Legendre functions and ( r) are determined 
from the radial differential equations 
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and 

G^( r=0) finite, G^(r=a) = 1 

The parameters u and e are the magnetic permeability and 
electric permittivity, respectively, and a is the electrical 
conductivity which is assumed to be an arbitrary function of 
radius. For the present, work , free space values for p and e 
have been assumed and MKS units are used throughout. 


For completeness we also present the spherical harmonic 
expansion for the TE portion of the incident magnetic field 
(aside from the factor e lu)t ) 

H r = H 0 sin, X '2fnr ^ +l > h (-T 1 ) P 1 < co58 > 
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(7) 


The transfer functions for the tangential magnetic field 
components can now be defined by decomposing the forcing field 
(1) into its components in the spherical polar coordinate system 
and forming the ratio, component by component, of the total 
tangential field components (3) to those of the forcing field. 
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The magnitude of the transfer functions for the 9 and cp 
components at the lunar surface (r=a) are 


{. 


cp 
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( 8 ) 


The transfer functions depend on the colatitude angle 9 but 
are independent of azimuth cp. We now define the point where 
9 = tt as the sub-k point, i.e., the point at which the wave 
vector is normal to and into the lunar surface. At the sub-k 
point 

^ dp|(cos9) p|(cos9) 

cos0 <T9 sin 9 ‘ ’ 

and therefore 

T 0 (9-rr) = ye-n) 

Sonett et al . [1971a] have defined a transfer function as 
the ratio, term by term, of the tangential magnetic field given 
by (3) to the driving field given by (7). This definition pro- 
vides a transfer function for each harmonic which is independent 
of 9 and cp and which is the same for both tangential components 
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It is difficult to interpret the experimental data in terms of 

9 

the transfer functions T because of the present inability to 
separate the data into individual harmonic contributions. We 
will see that for f < 0.01 Hz, « T Q « T 

~ 0 cp 


Approximations to the transfer functions T g and have 
appeared in the literature. Blank and Sill [1969] assumed a 
spatially uniform forcing magnetic field of the form H = mHge' 


i <ju t 


(in our notation) and thus considered a transfer function, which 
we refer to as Tg , that can be obtained from (8) by taking the 
limit 2iraA - 0 . The approximate transfer function Tg is re- 
lated to by 
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( 10 ) 


The approximation Tg is valid at low frequencies, but at 
frequencies greater than 0.01 Hz we will show that it departs 
significantly from (8). Since the experimental Apollo 12 LSM- 
Explorer 35 transfer functions are known at frequencies as 
high as 0.04 Hz, (8) is presently being used in the inversion 
of that data. 
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Three simple lunar electrical conductivity models will be 
used to illustrate the angular dependence of the tangential 
magnetic field transfer functions and to reveal the inadequacies 
of the approximation Tq at frequencies greater than 0.01 Hz. 

The first is a two layer model with a perfectly conducting core 
of radius b surrounded by an insulating shell. The second model 
is also a two layer moon but with a fi-nitely conducting core. 

The third is a moon of zero conductivity except for a thin con- 
ducting shell of conductivity o and thickness 6 centered at r=b. 
This model is suggested by the results of Sonett et al . [1971a, b], 
who inferred a near surface "spike" in electrical conductivity 
by inverting the combined Apollo 12 LSM and Explorer 35 harmonic 
data. To further simplify this model we assume that 5-0 while 
a6 remains constant; the thin conducting shell is thus represented 
by surface currents at r=b. 
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TWO LAYER MOON WITH A PERFECTLY CONDUCTING 
CORE AND INSULATING SHELL 


For this model, the contribution to the transfer function 
(8) of the radially dependent term is 



r=a 


1 X l±l + £(b/a) 2 ^ +1 ) 

a 1 ] . (b/a) 2l+] J 


( 11 ) 


If the conducting core is large, i.e. l-b/a«l, the right 
hand side of (11) reduces to (a-b) a result which is inde- 
pendent of l and thus allows the following numerical example 
to be scaled to other values of b. 


The transfer functions T Q and T have been calculated 

t) Cp 

for b = 1505 km over the frequency range 0 ^ f ^ 0.04 Hz. The 
upper frequency corresponds to the highest frequency for which 
transfer functions are obtainable in the Apollo 12 LSM-Explorer 
35 magnetometer experiment and the choice of b was motivated 
by the radial location of the near surface conductivity "spike" 
reported in Sonett et al . [1971a, b] . The computations assumed 
a wave propagation velocity of 400 km/sec and were carried out 
for a range of angles from 8 = tt, the sub-k point, to 0 = 2rr/3 
Calculations were not carried out for 0 < 2tt/3 since the 
assumption of induced field line confinement is expected to be 
best realized near the sub-solar point which we assume to be 
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near the sub-k point. Results for Tg are shown in Figure 3- 
For the values of 6 investigated here, is rather insensitive 
to changes in 0, i.e. (0=tt) = Tq(9=tt). The transfer 

functions Tg show a rolloff with increasing frequency for 
0 Si 130°, the rolloff being greatest at the sub-k point. The 
transfer functions do not approach unity as f-»0 because of the 
assumption of infinite core conductivity. For any finite value 
of electrical conductivity the transfer functions indeed approach 
unity as f-*0. Since Tg is relatively sensitive to 0, while T^ 
is not, one should expect a difference between the values of the 
two experimental transfer functions for LSM locations not coin- 
cident with the sub-k point. The transfer function T^ (also 
shown in Figure 3) is a monoton i ca 1 1 y increasing function. The 
approximate transfer function Tg is independent of both angle 
and frequency and is given by 



1 + \ (b/a) 3 
1 - (b/a) 3 


( 12 ) 


This is the low frequency limit of T^ . For the case illustrated 
in Figure 3, (12) would give a straight line at 3.75. 



TWO LAYER MOON WITH CORE OF FINITE CONDUCTIVITY a c 
AND AN INSULATING SHELL 


The radially dependent term for this model is 


(£) 


r=a 


a 


1 + 1 + l(b/a) 2l+] Bj, ( kb) 


1 - (b/a) 2 ^ +l B^(k c b) 


( 13 ) 


2 2 

where k = k (a=cr ), r = b is the core radius, and 
c v c J 


Bit kb) = 1 - 


(2t+l) j t ( k^b) 


k c b 4-1 (k c b) 


(Ik) 


Both T q and T were calculated as functions of f and 9 for the 
b cp 

case b = 1 566 km and c = 1 0 '' mhos/m, the best 2- layer fit to 

c 

the experimental transfer functions reported by Sonett et al . 
[1971a]. The results for Tg are shown in Figure 4. As in the 
previous case T^ showed only minor variation with 9 for tt:s9< 2tt/3- 
The function T 0 increases from unity at f = 0 and rolls off 
above 0.03 Hz for 9 > 1 40° . For 120° < 9 < 140°, T 0 is mono- 
tonically increasing. As in the previous case we expect a 
difference between the two experimental transfer functions when 
LSM is not at the sub-k point. The function T^ (the broad dashed 
line in Figure 4) is again monoton i ca 1 1 y increasing. The approxi- 
mation Tq , also shown in Figure 4, tends to lie in the region 
of the figure which corresponds to T 0 for 9 = 130° and is inaccurate 


at f ^ 0.01 Hz. 
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NONCONDUCTING MOON WITH THIN CURRENT LAYER 


Sonett et al . [1971a, b] have reported a near surface 
shell of high electrical conductivity imbedded within a 
relatively insulating Moon except perhaps for a small conducting 
core. Thus, transient eddy currents in the Moon are concentrated 
mainly in this highly conducting thin ‘shell. A mathematical 
'idealization of these findings is a nonconducting Moon contain- 
ing a current layer of negligible thickness situated at r = b and 
characterized by a nonzero value of the product crS . 

2 -£ jf. j . 1 

Solutions of (5) with k ^ 0 are proportional to r , r . 

-X 

In the region 0 ^ r < b the solution proportional to r must 
be discarded and we find 



From boundary condition (6) and the continuity of the radial 
component of the magnetic field at r = b we obtain 


C t a~ l + a l+] =1 , (16) 

\ b ' 2< "' + 


(17) 
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The additional boundary condition required to determine the 
solution follows from the fact that the jump in the tangential 
magnetic field is equal to the surface current density at r = b 

1 i m Hg ( r=b+ -^) - Hg(r=b- -^ ) = cr6 E ( r=b) . (18) 

6~ > 0 

For the TE mode, the electric field component E^ may be found 
by integrating the 0-component of the Maxwell equation vx_E = - 
3j3/3t and noting that E r = 0 and Hg is given by (3). The result 
i s 

= 0rr -, dp! ( COS0 ) 

E v = ' ^ fH 0 sintp L 7 h ( V> h — —55 ’ (,9) 

•1=1 

and from (3), (15), (18) and (19) we obtain 

U+1)D^ - A t {l+] - 3 «uut] , (20) 

where t = ap6b/3 > (21) 

has been shown to be the appropriate time constant (Cowling 
decay time) of the current layer model for response in vacuum 
to an abrupt change in a spatially uniform magnetic field 
[Schubert and Colburn, 19711. Combining ( 1 6 ) , (17), and (20) 
we find that the radially dependent term, evaluated at r=a , is 
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]_ 

a 


h ^ 

1 + ] + B^( 3 iujT ) 

, 21 + 1 

1 “ (-|) B ^( 3-i lut ) 


where B^(3iu)r) = (l - ^+1) 


( 22 ) 


(23) 


The radially dependent term defined by (22) is formally 
identical to the one for the two layer Moon with a conduct- 
ing core and an insulating shell, (13)- The difference between 

the two lies in the definitions of B.(3Iuut) and B ( k b). If 

"C c 

(jut and k c b become infinite, both (13) and (22) reduce to (11), 
the radially dependent term for the two layer Moon with a 
perfectly conducting core and an insulating shell. 

The illustrative model has the conducting layer at 
b = 1505 km and the product a 6 = 150 mhos. This latter quantity 
corresponds, for example, to a 50 km wide layer with a con- 
ductivity of 3x10 mhos/m. The approximation Tg and the trans- 
fer functions T Q and T^ for this case are shown in Figure 5. 

As in the previous two cases, T^ is relatively insensitive to 
variations in 8 and is not shown in the Figure. The transfer 
function T 0 is sensitive to 0 variations and also exhibits two 
characteristics which are significantly different from those 
of the two layer model with a finitely conducting core (Figure 
4). First, consider the low frequency (f £ 0.0075 Hz) region 
of the first curves. Both transfer functions (Figures 4 and 5) 
reach a magnitude of 3 at f = 0.00625 Hz. The current layer 
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response is almost linear between f = 0 and f = 0.00625 Hz 
whereas the two layer model shows a much more rapid rise which 
is turning over when Tg = 3. This difference was pointed out 
by Sonett et al . [1971a] in their discussion of the inadequacies 
of models other than their proposed conductivity profile. The 
second characteristic of the current layer model is the 
flattening of its response at a much lower frequency, f ~ 0.015 Hz, 
as compared to the two layer model which just barely flattens 
at f = 0.03 Hz. In Figure 5, 

t - (i + tJ)}* 

from Sonett et al. [1971a, b], is shown with its error bars. The 
simple current layer model provides an excellent fit to the data 
at the lowest frequencies and is well within the error bars over 
most of the range 0.001 ^ f ^ 0.04 Hz. It should be noted that 
the data points represent averages of observations at presently 
unknown values of 9 , since identification of 0 would require 
empirical knowledge of the distribution of the wave propagation 
vectors of the disturbances forcing the lunar induction. 
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CONCLUSIONS 

It has been shown that there are two distinct transfer 
functions corresponding to the two orthogonal magnet i c f i el d 
components tangent to the Moon's surface; they are identical 
at the sub-k point but exhibit different functional behavior 

on the colatitude angle 0. One transfer function, T / is 

cp 

relatively insensitive to 0 while the other, Tg/ shows marked 

change as the colatitude angle decreases. Both functions 

show a significant high frequency rolloff for 0 > tt but Tg 

r 

changes to a monoton i ca 1 1 y increasing function between 9 = tt 
2 

and 0 = -^ rr . The 1 imitations of the low frequency approxima- 
tion Tg have been made quantitative by the exact theory of this 
paper. 

The present calculations provide an explanation for one 
of the unexpected features of the experimental transfer functions, 
namely a rolloff at the high frequency (f > .025 Hz) end of the 
spectrum (see Figure 1). This rolloff could not be understood 
on the basis of the approximation Tg to the theoretical transfer 
functions. Neither does the transfer function T^ exhibit a 
high frequency rolloff. The analysis of this paper shows that 
Tg and T^ do roll off at high frequencies if the observer is 
near the sub-k point. This phenomenon can be understood as 
follows. The relatively large magnitudes of the lunar transfer 
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functions result from the compression of induced magnetic 
fields into a shell between the confining solar wind and an 
inner "core" which excludes field lines at sufficiently high 
frequency. All the harmonics in (8) are thus amplified, how- 
ever, the effect is of decreasing importance as one proceeds to 

o 

higher harmonics, i.e. i decreases with increasing t at fixed 
f. Since the higher harmonics contribute proportionately more 
to the field at higher frequencies, i.e. since the power in 
the field is redistributed so that at higher frequencies 
relatively more power resides in the higher harmonics, the net 
effect is a rolloff in the amplification or transfer function. 

That the rolloff is a function of proximity to the sub-k 
point is also clear from (8) if one notes that the contribution 

1 dP I 1 1 

of each harmonic is proportional to either cos Q TPT~ or 'sTnl) P -t 
For t = 1 these factors are both unity whereas for 1 = 3 (the 
harmonics l = 1 and 1=3 are approximately in phase) they are 
(A - 15 sin 2 0) and ( 4 - 5 sin 2 0), respectively. As 0 de- 
creases from the value n , both these factors decrease, the first 
more rapidly, which accounts for the sensitivity of Tq on 0 
and the relative insensitivity of on 0. 

One final point concerns the observed asymmetry in the 
experimental data. Sonett et al . [1971a] showed data in which 
the Z-transfer function, corresponding to local north or 
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approximatel y perpendicular to the ecliptic, was in general, 
greater than the Y-transfer function, local east (Figure 1). 

The Z-transfer function corresponds to in this paper and the 
Y-transfer function corresponds to Tg (assuming the wave vector 
lies in the equatorial plane). From the theoretical analysis, 
based only on a TE mode interaction, Tg £ T^ and at frequencies 
below 0.01 Hz, Tg « T^ for all values ‘of 0. The asymmetry in 
the data is nevertheless present even at the lowest frequencies, 
which indicates that the asymmetry must be attributed to some 
other cause, for example the presence of a TM mode or an 
asymmetry in the lunar conductivity distribution. Because of 
the way in which Tg depends on f and 0 we expect the basic 
asymmetry to decrease with frequency, a result confirmed by 
the data ( F i gure 1 ) . 
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FIGURE CAPTIONS 


Figure 1 . 


Figure 2. 


Figure 3 . 


Figure A . 


Magnitude of the Apollo 12-Explorer 35 tangential 
transfer functions vs. frequency (from Sonett et al., 
1971a). The Y component is east and the Z component 
north at the Apollo 12 site. Note the rolloff at 
the high frequencies and the differences in the trans- 
fer functions for the two tangential components. 

The geometry of the lunar induction problem. The 
wave propagation direction of the external disturbance 
forcing the lunar response is in the C-direction. 

The angle 6 is the colatitude of LSM with respect to 
the direction of the wave vector. 

The transfer function Tg and the approximation T^ vs. 
frequency for a model Moon with an insulating shell 
and a core of infinite conductivity with radius 1505 km. 

The transfer functions Tg and T^ and the approximation 
Tq vs. frequency for a model Moon with an insulating 
shell and a core of radius 1566 km and uniform con- 
ductivity 10 mhos/m. Note the rolloff in Tg at 
high frequency for 6' > 1 40° . The approximation Tg 
is inaccurate at f s 0.01 Hz. 



gure 5. 


Theoretical transfer functions for a model Moon 
which is insulating except for a current layer at 
1505 km with u6 = 150 mhos. The data shown (points 
with associated standard deviations) are values of 

t ■ a < t z + T ?>r 

from Sonett et al . [1971a, b] . The current layer 
model provides an excellent fit to the data at the 
lowest frequencies and a reasonable fit at the high 
frequencies considering that the data represent 
averages of observations at unknown values of 8 . 
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